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A simple model to link the properties of quasars to the 
properties of dark matter halos out to high redshift 
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ABSTRACT 

We present a simple model of how quasars occupy dark matter halos from z = to 
z = 5 using the observed meH — c relation and quasar luminosity functions. This pro- 
vides a way for observers to statistically infer host halo masses for quasar observations 
using luminosity and redshift alone. Our model is deliberately simple and sidesteps 
any need to explicitly describe the physics. In spite of its simplicity, the model re- 
produces many key observations and has predictive power: 1) model quasars have the 
correct luminosity function (by construction) and spatial clustering (by consequence); 
2) we predict high redshift quasars of a given luminosity live in less massive dark 
matter halos than the same luminosity quasars at low redshifts; 3) we predict a factor 
of ~ 5 more IQ^-^Mq black holes at z ~ 2 than is currently observed; 4) we predict a 
factor of ~ 20 evolution in the amplitude of the men — -^haio relation between z = 5 
and the present day; 5) we expect luminosity dependent quasar lifetimes of between 
tq ^ 10^~® yr, but which may become as short as 10^~^ yr for quasars brighter than 
L*; 6) while little luminosity dependent clustering evolution is expected at z 1, in- 
creasingly strong evolution is predicted for L > L* quasars at higher redshifts. These 
last two results arise from the narrowing distribution of halo masses that quasars 
occupy as the Universe ages. We also deconstruct both "downsizing" and "upsizing" 
trends predicted by the model at different redshifts and space densities. Importantly, 
this work illustrates how current observations cannot distinguish between more compli- 
cated physically motivated quasar models and our simple phenomenological approach. 
It highlights the opportunities such methodologies provide. 
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1 INTRODUCTION 

Quasars represent a unique population of objects in the 
Universe that encapsulate many otherwise diverse areas 
of physics. These include extreme environments of gravity 
(black holes), sub-to-kiloparsec-scale dynamics (black hole 
two and three body interactions and host galaxy merg- 
ers or secular triggers of quasar activity) , sub-to-kiloparsec- 
scale hydrodynamics (gas infall, accretion disks and quasar 
winds), and quasars as cosmological probes of the evolving 
large-scale cosmic web. They are among the most luminous 
objects in the Universe. The energy liberated during a single 
quasar event can outshine the entire stellar light of the host 
galaxy. After fading, their presence can still be measured 
through the local quiescent black hole population. 

Although much work has been done to describe the 
physics of black holes and their evolution, the majority of 
what we know remains primarily phenomenological. Two 
fundamental correlations are observed: the ttibh — o" relation 



Tibuigc relation ('Ferrarese & Merritt! '2000|: 
Tremaine et al. 2002; Marconi & Hunt 2003; Haring & Rix 



and the msH 



20041 ). Although initially surprising (why should the sub- 



parsec physics of black hole growth correlate with the kilo- 
parsec properties of the galactic bulge?), it is now be- 
lieved that these relationships simply reflect t he physics 
of a common form ation mechanism (ISilk fc Rees 199^; 
iHopkins et al.|[2006ah . For example, galaxy major mergers 
may simultaneously drive growth in the bulge and force gas 
into the central regions of the galaxy to fuel the black hole 
(and hence a quasar). Secular processes may be operating 
to sim ilar effect, such as bar instabilities (|Sellwood fc Moord 
Il999l ). 

There is a need to understand the phenomenology of 
black holes and galaxies in greater detail. Active black holes 
are thought to have great impact on the evolution of their 
hosts. During a quasar event, rapid hole growth occurs and 
outflows drive winds that may liberate gas f r om the galaxy 
(|Di Matteo et al.ll2005l : iHopkins et"aLll2006al : iThacker etall 



2 D. J. Croton 



I2OO6I ). Such gas is considered fuel for future star formation, 
and without it new stars no longer form and the host galaxy 
subsequently reddens and fades, a kind of galactic extinc- 
tion. However, supermassive black holes do not always die 
with their galaxy, but are often later found in a low luminos- 
ity state. Heating from low luminosity active galactic nuclei 
(AGN) provide a long term energy source that suspends 
the cooling of gas from the surrounding hot halo (the so 
called "radio mode" solution to the " cooling flow" problem: 
ICroton et alll2006l : [Bower et aLlbOOd ). This heating/cooling 
balance is known to maintain the aged appearance of many 
massive local galaxies. Understanding AGN and black holes 
has now become essential to understanding galaxies, and 
hence modelling the co-evolution of both has become a sub- 
ject of great interest. 

In this paper we build a general model of how black 
holes and quasars occupy dark matter halos and how this 
occupatio n evolves with time. Our met hod i s similar in spirit 
to tha t bv iMarinoni fc HudsonI (|2002l ) and IVale fc Ostritol 
l|2004l ) but for quasars rather than galaxies. We work from 
a minimal set of assumptions and use two key observational 
constraints: the man — o" relation and the quasar luminosity 
function. Under such constraints the model naturally repro- 
duces many quasar and black hole observations out to red- 
shifts as distant as z ^ 5. A number of predictions are given. 
This model is deliberately simple; it sidesteps any attempt 
to explicitly describe the quasar triggering mechanism, the 
details of black hole accretion, or the hydrodynamics of sub- 
sequent quasar winds and outflows. Although understanding 
such detail is certainly desirable, we show that the current 
observations do not discriminate between the more com- 
plicated physically motivated modelling of quasars and our 
simple phenomenological model. 

The outline of this paper is as follows. In Section [5] we 
describe the construction of the model and the observational 
data used to constrain it. In Section [3] we explore the vari- 
ous consequences of this model, comparing to observations 
where available and making predictions where not. Section |4] 
provides some discussion, placing the model into a broader 
context of black hole and galaxy co-evolution. Finally, Sec- 
tion [5] summarises our main results. Unless otherwise stated, 
we assume a standard WMAP first year ACDM cosmology 
with SIa — 0.75, Qm ~ 0.2 5, as = 0.9 and Hu b ble constant 
Hp = 100/ikm/sMpc-i jSpergel et all |2003| ; [Seiiak et all 
I2OO5I ). We choose the value of the Hubble parameter, h, to 
be either h — 0.7 or ft = 1.0 depending on the context; this 
will be clearly marked. 



2 ACCURATELY POPULATING DARK 
MATTER HALOS WITH QUASARS 

We build our quasar model in two parts. First, we map 
quasar luminosity onto dark matter halos using the halo 
virial properties. Second, we determine which of these ha- 
los actually host a quasar as a function of luminosity at 
any given redshift. Both parts are undertaken using obser- 
vational constraints only, notably the uibh — a relation and 
quasar luminosity function. 

To begin we require knowledge of the dark matter halo 
population and its evolution. There are several ways this 
can be achieved at a given redshift for a given cosmol- 



ogy. Analytic me t hods, such as those rec e ntly d escribed in 
iNeistein fc Dekell HqoI) and lZhang et al.l (120081 ). are not as 
useful here as they produce halo merger trees lacking spatial 
and velocity information. We will later need these properties 
in our analysis. 

Instead, we turn to a numerical N-body simula- 
tion of dark matte r evolu tion, the Millennium Simula- 
tion ISpringel et al.1 l2005d ). This simulation, run using 
a WMAP1-I-2DFGRS cosmology, follows the evolution of 
10 billion dark matter particles in a box of side-length 
500 Mpc from 2: = 127 to 2: = 0. Within the simulation, 
both halos and subhalos (i.e. the bound sub-structure within 
a given halo) are accurately resolved down to virial masses 
of less than lO" Mq, more than sufficient for our purposes 
here. Note that in what follows we do not discriminate be- 
tween halos and subhalos when populating the simulation 
with quasar^. Thus, it is possible in our model for given 
halo to host more than one quasar at any given time. Small- 
scale clustering measures of quasars indicate that this may 
indeed be the case (e.g. Hennawi et al.l I2OO6I : iMvers et al.l 
l2008l : iPadmanabhan et al.ll2008h . although it is not critical 
for our current analysis. 



2.1 Linking quasar luminosity with halo mass 

We start with the dark matter halo virial mass, Mvir, of 
each Millennium Simulation halo. Our goal is to first relate 
its properties to the expected velocity dispersion of an occu- 
pying galaxy, a, and then to the quasar luminosity through 
the rriBH — o" correlation. Readers who are only interested in 
the final relations should skip ahead to the equation sum- 
mary in Section 12.1.11 (see also equation [13] in Section 14. 2p . 

Dark matter halos are identified at each redshift as re- 
gions of the simulation whose mean density inside a spherical 
aperture exceeds 200 times the critical density of the Uni- 
verse. The virial mass of a halo and its virial velocity, Vvir, 
are then related by 

Kir= [10GJi'(z)A/virl'/' , (1) 

where Mvir has units of H'^Mq, G — 4.3 x 
10~^ (km/s)^ Mpc Mq ^ is Newtons gravitational constant, 
and 

H{z) = HoE{z) = Ho [an(l + zf + f^A]'/' (2) 

is the value of the Hubble constant at redshift z (lHogelll999l '). 
Here, Ho = 100 /ikm/sMpc~^ is the local value of the Hub- 
ble constant (with h the dimensionless Hubble value), while 
f2m and Qa are the universal mass and dark energy densi- 
ties. Note, that both A/vir and H{z) above must share the 
same h value for dimensional consistency. 

Halo virial velocity is typically related to the observed 
galaxy circular velocity, Vc, by 



(3) 



with 7 a parameter of order unity (see Section [2. 3|l . This in 
turn can be related to the velocity d ispersion of the g alaxy, 
a, through the observed correlation (|Baes et al.ll2003 



^ Throughout this paper our use of the term "halo" includes 
both halos (i.e. quasars within central galaxies) and subhalos (i.e. 
quasars within satellite galaxies). 



Linking quasars to dark matter halos 3 



log (— ) = (-0.22 ± 0.03) + (1.04 ± 0.12) log ( — ) , (4) 

where uq — 200km/s. Equations [l]|4] connect the virial mass 
of a dark matter halo at redshift z to the velocity dispersion 
of the occupying galaxy. 

A well defined correlation between stellar velocity dis- 
persion and black hole mass is observed in the local Universe, 
the mBH — o" relation (jTremaine et al.ll2002l ): 



log I 



(8.13 ± 0.09) + (4.02 ± 0.44) log (— ) , (5) 



where /170 is the Hubble parameter with Ho = 
70 km/s Mpc~^, and ag = 200km/s. For practical purposes, 
when using equation [5] in our model we include the observed 
0.3 dex dispersion in mBH (for simplicity we apply the same 
dispersion at all values of a). We assume no evolution with 
redshift in either amplitude or slope of the relation. We dis- 
cuss this further in Section [3.21 

With a value of men for each dark matter halo we can 
now determine the luminosity of a quasar that may occupy 
the halo. We take this (bolometric) luminosity as some frac- 
tion, r;, of the Eddington luminosity: 

Lq = ri Lcdd 

= 7? 3.3 X 10* mBH {h^o Lq) 

= 7? 1.3 X 10^* mBH {hro erg/s) . (6) 

Equation |6] completes our sought after connection between 
virial mass and quasar luminosity, which we will now sum- 
marise. 



2.1.1 A summary of the important equations and 
relationships 

Equations [T]|6] provide a mapping between halo virial mass 
and quasar bolometric luminosity. These equations can be 
reduced to the following relation 



log (^^^) 
\ 1012 hz} Lo / 



-1.99 ±0.33) 



+ (1.39 ±0.22) log 



Hiz)(^ 



Mvir 



10" /i-i M, 



(7) 



where H{z) is defined by equation[2] Note that, for internal 
consistency, H{z) must be calculated using the same value 
of h as Mvir, but for Lq we have assumed the standard ob- 
server value of ft = 0.7. Quasar luminosity here is dependent 
on both parameters 7, describing the relationship between 
virial and circular velocities (equation[3ll , and rj, the Edding- 
ton luminosity fraction that quasars are assumed to shine at 
(equation [B]). Observational errors from equations |3] and [S] 
have been propagated throughout; their impact on the re- 
sults are discussed in Section [4.31 

Equations IH5I can also be reduced to map directly be- 
tween halo mass and black hole mass 

^ (-2.66 ±0.33) 



lo: 



108 /i-i Mq 



+ (1.39 ±0.22) log 



lH{z) 



Mvir 



1013 /i-i Me 



(8) 



This equation contains a redshift dependence through H{z), 
something we will investigate later in Section [3.61 

Finally, we provide a few useful equations from the lit- 
erature for converting to different optical filters. First, to 



convert between bol ometric quasar lumi nosity and 6j-band 
absolute magnitude, [Croom et al.l l|2005h provide 



-2.66 log(LQ) ±79.42 



(9) 



where Lq here is in watts and assumes h = 0.7. Second, 
to convert between B-band, 6j-band (used by the 2dF QSO 
Redshift Survey, hereafter 2QZ), and i-band (the Sloan Dig- 
ital Sky Survey Quasar Survey, hereafter SDSS) (k-corrected 
to 2 = 2) quasar lum inosit ies we use the convers ions given 
in lCroom etlo] l|2005h and lRichards et al.1 l|2006D 



Mb = Mftj ± 0.06 , 
Mi[z = 2] = Mb, -0.71 . 

2.2 Deciding which halos host quasars 



(10) 



Although Section [2.1 [ provides a mapping between halo mass 
and quasar luminosity, not all halos will actually host a 
quasar at any given time. To decide which do we use the 
quasar luminosity function. 

Specifically, we constrain the model to have the correct 
luminosity function at a giv en epoch usin g the technique 
of abundance matching (e.g. IConrov et al.l 2006 and refer- 
ences therein). First, we calculate the observed cumulative 
luminosity function, from bright to faint quasars, which pro- 
vides a smoother representation of the data than the lumi- 
nosity function alone. Then, starting at the bright-end, we 
move faint-ward down the cumulative luminosity function 
in magnitude bins of AM = 0.01 and randomly dilute the 
number density of model quasars in the same magnitude bin 
to ensure they match the same cumulative abundance as the 
observations. Hence, by construction, the model reproduces 
the observed quasar luminosity function. 

The observed quasar luminosity function is only mea- 
sured at discrete redshifts, and we would like to be able to 
build our model at any arbitrary redshift. Hence, instead of 
using the data itself, we abundance match to a functional 
fit of the data across multiple redshift s. The two da t a sets 
we will later compare with are those of Croom et al.l ( 2004h 
(2QZ), covering 0.4 < 2 < 2.1, and iRichards et al.l (|2006h 
(SP SS), covering 0.5 < z < 4.8. 

ICroom et al.l (|2004l ) model their data using an evolving 
double power-law with the form 



$(Lb,,z) 



(L,,/L*J-±(L,,/L-^)-/3 ' 
ic lumino 

Ll,{z)^Ll,{z = 0) 10*1 -+'=2-' 



(11) 



where the characteristic luminosity, L^^, is a function of red- 
shift 



(12) 



The authors used this functional form to fit the 2QZ quasar 

luminosity function out to z ~ 2^ 

Interestingly however, the ICroom et al.l double power- 
law is a good fit to the quasar luminosity function well be- 
y ond its original fi tting range. This can be seen i n figure 19 
of lRichards et al.1 l|200d ). where the|c room et al.l result con- 
tinues to provide a good match to the SDSS quasar lumi- 
nosity function out to redshift threfl Hence, to maintain 



2 Note that [Richards' et al.l fit their SDSS data using a simpler 
single power-law, which does not realistically represent the quasar 
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Table 1. Parameters for the do uble power-law fit to the joint 2QZ and SDSS quasar luminosity functions at z < 3 and z 3 fSection l2.2l 
and eauations llll and ll2l l. As in lCroom et alj ll2004f) , h = 0.7 was assumed in the fit where relevant. 



redshift range 


a (bright slope) 


/3 (faint slope) 


M*^iz = 0) 


fei 


fc2 


*,(Mpc~^mag~l) 


z < 3.0 


-3.31 


-1.09 


-21.61 


1.39 


-0.29 


1.67 X 10-^ 


z ^ 3.0 


-3.31 + 0.5(2 - 3) 


-1.09 


-21.61 


1.22 


-0.23 


1.67 X 10-" 



the most realistic shape of the quasar luminosity function 
across the widest p ossible redshift range we adopt the fit of 
ICroom et al.l (|2004l ) to z = 3. Above this redshift we mod- 
ify their parameters slightly, adding evolution to bright-end 
power- law slope and softening the decline of Ll_^ with red- 
shift. The fitting values we adopt are given in Table [1] The 
fits are shown below in Section [3.31 

2.3 Model assumptions and parameters 

In order to keep our model as simple as possible while re- 
maining accurate and (most importantly) understandable 
we list below our main simplifying assumptions: 

• Both the local rriBH — o" relation (equation [5} and the 
local correlation between a and Vc (equation [Jjl exist at all 
redshifts and do not evolve with time. 

• The relationship between halo virial velocity and galaxy 
circular velocity can be described by equat ion [3] with 7 — 
1.0 (see Appendix A5 of lPorciani et al.ll2004l for a discussion 
on this point). 

• Quasars can be modelled as a simple "light bulb", 
where, at any given time, they are either on or off. 

• "On" quasars shine at half the Eddington luminosity 
(i.e. 7] = 0.5 in equation[6l). This is the mean observed value 
measured bv lMcLure fc Dunlod ||2004) during the hei ght of 
quasar activity at z = 2 (see also iMarconi et al]|2004) . 

The impact of the uncertainty of these assumptions on our 
results is further discussed in Section [4.31 

Although the above choices are reasonable they are not 
necessarily correct in detail. For example, quasar luminosi- 
ties are not either "on" or "off" , but follow a hght curve with 
a peak luminosity that is likely dependent on the specifics 
of the quasa r trigger, t he properties of the host galaxy, and 
redshift (e.g. lHopkins e t al. 2006aj). More complex modelling 
is required to capture these processes. Our goal here is not 
to model the specifics of the quasar population in detail, 
but rather we aim to find the simplest model possible that 
can match a set of key observations, in spite of the missing 
detail. We will discuss this further in Section [l] 



3 CONSTRAINTS, CONSEQUENCES AND 
PREDICTIONS TO HIGH REDSHIFT 

3.1 The quasar luminosity- halo mass relation 

The relationship between quasar luminosity and dark matter 
virial mass is given by equation [T] This is plotted in figure [T] 

faint-end. T heir data does not pro be below L* at these redshifts 
(see however iFontanot et al.ll2007l) . 




12.0 12.5 13.0 13.5 14.0 

log Mhoio (h^'Mo) 

Figure 1. The relationship between quasar luminosity and dark 
halo virial mass at various redshifts from 2 = to z = 5 (equa- 
tionO. Magnitudes are provided in both the bj-band used by the 
2QZ survey and i-band used by the SDSS (k-corrected to 2 = 2 - 
see equation IIOI I . Note that quasars of a fixed halo mass become 
brighter with increasing look-back time. 



for select redshifts out to z — 5. On the left axis we plot 
&j-band quasar absolute magnitude for comparison with the 
2QZ survey results, and on the right axis we show i-band 
magnitudes for comparison with the SDSS quasar results 
(see equation [10} . 

Figure [T] reveals that model quasars hosted by halos of 
a given mass get brighter the further back in time you look. 
For example, in halos of lO^^'^M©, quasars brighten in lumi- 
nosity by about two magnitudes between z = and z = 2, 
and by another one-and-a-half magnitudes between z = 2 
and z = 5. One can see this by eye directly in the observa- 
tional data, ho wever, without deferr ing to theory or models. 
The results of ICroom et al.l (|2004l ) show similar amounts 
of brightening in the 2QZ quasar luminosity function be- 
tween z = . 4 and z — 2.1, while for the same catalogue, 
ICroom et all (l2005l ) use clustering to infer that the masses 
of L« quasar dark halo hosts across the same redshift inter- 
val remain approximately constant at IO'^^-'^^Mq. In this 
sense, our model is simply mimicking the data from which 
it was constrained. 

From a theoretical point-of-view, the evolution in am- 
plitude seen in figure [T] originates from a redshift depen- 
dence in equation [TJ where, at fixed halo mass, the virial 
velocity of a halo incr eases with increasing redshift (see also 
IWvithe fc Loebll2003l ). This behaviour simply tells us that 
halos of equivalent mass live in higher cr-peaks at higher 
redshifts. 
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100 200 500 

dispersion (km/s) 

Figure 2. The rviBH — f relation between quasar black hole mass 
and galaxy velocity dispersion. We assume in the model that this 
relationship is redshift independent (s ee Section 12.1 [I. A fit to 
the observed local relation measured bv lTremaine et al IS 
shown by the solid thick line. 




-22 -24 -26 -28 

f^bj - 5logh,„ 



Figure 3. 2QZ quasar luminosity functions llCroom et al.ll2004l 
covering z = 0.40 to z = 2.10 (symbols with error-bars). The 
(diflicult to see) dotted lines show the fit to the data described 
by eQuation I 1 1 I with double power-law parameters given in table[T] 
Solid lines show the model luminosity functions determined from 
abundance matching to the observations (Section 12.21 1. 



3.2 The men — (y relation 

The man — cr relation is the primary observation used in 
our model to map quasar luminosity onto mass. To produce 
a realistic quasar mock catal ogue we include t he ob served 
0.3 dex of scatter reported in iTremaine et al] (|2002l ) when 
performing this mapping. Our model relation is plotted in 
figure [21 where the points are a randomly selected sample 
of quasars from the model, and the thick solid line shows 
the best fit to the observed correlation. It is important to 
note that, being a model quasar population, we are 100% 
complete down to low dispersion, unlike the measured data. 
Such measures, especially at low dispersion, still remain ob- 
servationally formidable. 

When constructing the model it was interesting to find 
that including the above dispersion does not appear im- 
portant for its success in any way. All results presented in 
this paper are essentially unchanged if we had rather sim- 
ply made a direct mapping between ttibh and a, ignoring 
the scatter. Evidence is emerging to suggest that a tight 
corre lation exists between Lq and Mhaio (e.g. IWhite et al.l 
I2OO7I ). Understanding this result will be critical when using 
quasar clustering measures to observationally constrain the 
host masses of high redshift quasars. 

We have made a critical assumption when constructing 
our model that the rriBH — c" does not evolve with redshift, ei- 
ther in amplitude or slope. An evolving relation would place 
quasars of a given magnitude in either more massive halos 
(for a decreasing amplitude with increasing redshift) or less 
massive halos (for an increasing amplitude with increasing 
redshift) . This change will be reflected in the clustering prop- 
erties of the quasars themselves. As we show in Section [3.71 
our non-evolving tubh — o" relation assumption produces a 
quasar population whose 2-point function matches the ob- 
servations extremely well out to at least z ~ 4. The obser- 
vational picture as to whether the meH — cr relation evolves 
remains unclear (see ICrotonI |2006| for a discussion) , so for 
the current work we retain the non-evolving assumption. 



3.3 The quasar luminosity function to z ~ 5 

As discussed in Section [2.21 we use abundance mapping to 
select those halos with quasars that will produce a luminos- 
ity function identical to that observed at each redshift. To 
maximise its versatility we match the model to a fit of the 
data that varies smoothly with redshift (Section l2.2l and ta- 
ble[T|, rather than to the data itself (which limits us to work 
at the observed redshifts only) . 

Model luminosity functions are shown for z < 2 
in figure [5] comparing with the 2QZ survey results of 
[Croom et all liooi), and fo r z > 2 in figure [H compar- 
ing with the SDSS results of [Richards et ahl (|2006l ). In each 
panel of both figures, the symbols show the observed data, 
the dotted line the double power-law flt to the data, and the 
solid line is the model. 

It is no coincidence that the model line and double 
power-law fits are essentially indistinguishable for all but 
the lowest space densities at the highest redshifts (at which 
we are limited by the finite size of th e Millennium Simu- 
lation). Note that our extension to the! Groom e t al. ] (12004 ) 
double power-law representation of the 2QZ luminosity func- 
tion (table [T]) provides a good fit to the data up to the high- 
est redshifts probed by the SDSS quasar survey in figure [4] 

3.4 Quasar lifetimes 

The quasar lifetime, tg, is defined here as Nq = tq/tH, 
where Nq is the fraction of halos at a given redshift that 
host a quasar, and tn is the Hubble time at that redshift. 
Nq is a function of both mass and redshift and is determined 
in our model through abundance matching model quasar lu- 
minosities to the observed quasar luminosity function (Sec- 
tion l2.2|l . It is, in this context, a normalised quasar selection 
function. 

In figure [SI we present quasar lifetimes as a function of 
limiting faint quasar magnitude for five redshifts, ranging 
from z = 0.5 to z — 4. Lifetimes in general range from be- 
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Figure 4. Similar to figure [S] but this ti me for quasars in th e 
SDSS survey covering z = 2.01 to z = 4.75 llRichards et aLl2006l) . 
Again, dotted fines indicate the double power-law fit to the data 
(equation 1111 and table [TJ, while the solid lines show the model 
quasar luminosity function (Section [2?2j. 

tween 10^ to 10* years. At lower redshifts we see a distinct 
turn-over in tg, where brighter quasars have much shorter 
lifetimes (as short as 10^~® years) than fainter quasars (lev- 
elling out at ~ lO'^'* years). This turn-over occurs at approx- 
imately L* in the quasar luminosity function (equation [12]). 
The turn-over is less pronounced for higher redshift quasars, 
flattening somewhat and even increasing to high luminosi- 
ties ai z — 4. 

Due to the simplicity and transparency of the model we 
know exactly why model quasar lifetimes behave in this way: 
quasars occupy a narrower range of halo masses at late times 
relative to early times (skip forward to figure [TD]to see this). 
Because of this, low redshift bright quasars are rare amongst 
the abundant ~ 10^^ Mq mass halos they occupy (hence the 
turnover at bright luminosities) , whereas high redshift bright 
quasars are frequent among their rare ~ mass halo 

hosts (hence tg remains flat). Fainter quasars (those with 
~ Lg < L*) tend to always commonly populate mostly 
abundant halos, again resulting in a relatively constant tg. 
We will return to this point in Section [4.21 

Physically speaking, it is important to realise that the 
trends seen in figure [S] do not result from the explicit mod- 
elling of a changing Ed dington accretion frac tion, as has 
often been explored (e.g. [Hopkins et "al]|2005bh . Quasars in 
our model are assumed to always accrete at a fixed fraction 
of the Eddington rate across the quasar hfetimejf]. This is a 

Note that when tg becomes longer than the doubling time for 
the black hole ( ^ 10* yr) black hole mass and quasar luminosity 
change non-trivially across the quasar lifetime. This is an inter- 
nal inconsistency that all light bulb models must navigate. For 
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Figure 5. Predictions for quasar lifetimes, tg, as a function of 
limiting faint quasar magnitude, at five different redshifts from 
z = 0.5 to z = 4.0. At 2 2 and Lg > L* the quasar lifetime 
shortens considerably. At higher redshifts, tg remains approxi- 
mately constant or may even increase for the brightest. 

key difference between our model and many previous works, 
and may explain its ability to simultaneo usly match such 
a wide range of observations. For example, IWvithe fc Loebl 
find an over-production of bright quasars at low red- 
shift. In our model, such quasars have very short lifetimes, 
and hence are not commonly seen in surveys. This may sim- 
ply arise due to the dwindling supply of cold gas in massive 
systems at late times |Fabianiil994li' ). 

3.5 The active and passive black hole mass 
functions 

Observations suggest that a black hole gains the majority 
of its mass while in th e active high accretion (quasar) phase 
(jHeckman et al.ll2004[ ). The cumulative effect of such mass 
growth over cosmic time is measurable in the local (passive) 
black hole population. Our model makes a prediction for 
the active black hole mass function. It is important to note 
that this prediction arises via the dual constraint of linking 
quasar luminosity to halo virial mass through the tubu — cr 
relation (equation [TJ, and from matching the abundance of 
black holes to the quasar luminosity function (Section l2.2|l at 
various redshifts. We did not tune the model in this regard 
to force a particular outcome. 

In figure [6] we show the observed cumulative black 
hole mass function for both z — passive and 2 = 2 ac- 
tive black holes. This fi gure is adopted from figure 6 of 
[McLure fc Dunlod i 20041 ) . The upper thin solid and dashed 
lines show the local observed results inferred from the 
TiBH— bulge luminosity relation and men — o" relation, re- 
spectively. The three data points show the cumulative SDSS 
quasar mass function at z = 2 for three different limiting 
black hole masses. The upward pointing arrows on each in- 
dicate that the space density measured in each bin is incom- 
plete, and hence provide only a lower limit to the true mass 
density. 

The three thick lines in figure [6] show the model pre- 
diction for our complete sample of quasar black holes at 

our work, this is most relevant for extreme luminosity quasars at 
the highest redshifts probed. Ignoring this effect for the sake of 
simplicity does not change our conclusions in any qualitative way. 
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Figure 6. The cumulative black hole mass function. The lower 
solid, dashed and dotted lines show the model result at three 
different redshifts, while the filled circles indicate lower lim- 
its measured from the S PSS z = 2 quasar mass function of 
iMcLure fc Dunlopl l l200'tf) . Note that the model predicts a fac- 
tor of ~ 5 more lO*'^ Mq black holes (quasars) than currently 
be seen in the data at 2 = 2. The upper solid and dashed lines 
show the local (passive) cumulative black hole mass function, cal- 
culated in two different ways. The passive function is built from 
the continual production of quasars across cosmic time. 



redshifts 0.5, 2.0, and 4.0, as indicated in the legend. For 
black holes more massive than log uibh ~ 9.0 at z = 2 our 
model is close to the lower limit found in the SDSS data. 
If true, the model indicates (perhaps unsurprisingly) that 
much of the massive end of the black hole mass function 
forms solely from accretion during this time of peak activ- 
ity. At lower black hole masses, logmBH ~ 8.5, our quasar 
model predicts an excess of fainter quasars yet to be seen in 
the z ~ 2 data. At higher redshift, z = 4, the model shows 
the continuing build-up of massive black holes but with less 
activity at lower masses. By late times, z = 0.5, massive 
black hole mass growth has largely stopped while the low 
mass holes continue to grow. 

This last behaviour is a black hole manifestation of 
the popular "downsizi ng" paradigm l|Heckman et al.ll2004l : 
iMerloni fc Heinj |2008l). It is interesting that such shift in 
black hole mass growth with time arises naturally from the 
model constraints alone. If we skip ahead to figure [TOl (Sec- 
tion l4.2p we can see the reason why. Here, dashed lines show 
the changing space density of halos/quasars with time, while 
horizontal dotted lines show black hole mass. At high red- 
shift, very low space density contours are mostly flat (or only 
slowly rising) and correspond to a fixed black hole mass of 
~ ]^g9-io^g lower redshifts all density contours turn 
over sharply, and black hole mass decreases by up to a few 
orders-of-magnitude at fixed space density. Hence, downsiz- 
ing is predicted for all objects ai z ^ 2. At z ^ 2 the model 
predicts that low space density objects (i.e. the most mas- 
sive) should show no downsizing trends (relative to 2: = 2), 
while higher space density objects (i.e. closer to L*) should 
be "upsizing" , especially above redshifts 2; ~ 4. 
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Figure 7. Evolution in the rrtBH — ^^halo relation and halo mass- 
to-bolometric luminosity ratio out to z = 5 (equation [Sj . There 
are two model predictions here: 1) that, at a give redshift, more 
massive halos should show an increasingly larger ttibh — ^hslo 
ratio (driven by accelerated black hole growth), and 2) the overall 
amplitude of this ratio should evolve with look-back time (see 
Section 13.611 . 



3.6 Evolution in the niBH Mhaio relation and 
mass-to-light ratios 

Section 12.11 and equation |8] provide an analytic prediction 
for how black hole and dark matter halo mass are related. 
This relationship includes a redshift dependence through 
E{z) = [ri,n(H- 2:)^ -t-SlA]^^^, implying that the niBH - Afvir 
ratio should evolve with time, with more massive black holes 
occupying dark matter halos of a fixed mass at higher red- 
shifts relative to lower redshift. 

In figure [7] we use equation |8] to plot the hole-to-halo 
mass relation at five different epochs, from 2 = to z = 5. 
At any given redshift the ratio increases with increasing 
halo mass, implying that black holes become proportionally 
larger the more massive the halo is. The evolution in the am- 
plitude of mBu/Mvir with redshift can also be clearly seen, 
with the ratio changing by a factor of ~ 5 between 2 = 
and 2 = 2, increasing to a factor of 20 by 2 = 5. Previous au- 
thors have attempted to quantify the chang e in black hole to 
host galaxy /halo properti e s with time (e .g. [Robertson et al.l 
I2OO6I : iMcLure et all 1200^ : ICrotonll2006l ). Evolution of this 
type is a clear predi ction of our model. Our results ar e simi- 
lar to those found bv lWvithe fc Padmanabhan ll2006ah using 
different techniques (see also IWvithe fc Loebl(200 3h. 

We can equivalently recast figure [7| as a changing mass- 
to-light ratio using Equation [T] The right axis in figure [7] 
shows this result. Halos with masses greater than 10^*Mq 
host quasars with mass-to-light ratios less than unity re- 
gardless of the redshift of interest. This is also true of all 
lower mass (Mv 
shifts 2^2. 



10 Mq) quasar/halo systems at red- 



3.7 Quasar clustering to z ~ 4 

The clustering of a given population of quasars will depend 
both on the masses of quasar hosts and the luminosity range 
which defines the quasar sample. Both key observations used 
to constrain our quasar model, the rriBH — o" relation and 
quasar luminosity function, are relevant to shape the model 
2-point function. 
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Figure 8. Model comparisons to tlie 2QZ (jPorciani et al.l[2004h and SDSS (|Slien et al.ll200'^ projected redsliift-space quasar correlation 
functions, from 2 ~ 1 to z 4. At all redshifts our model is a good match to the data (however with some small under-prediction of the 
clustering amplitude at z > 3). We also find a hint of steepening in the clustering on small scale scales (rp ^ l/i~^Mpc). 



In figure [8] we present the observed and model pro- 
jected correlation functions at five discrete redshifts from 
2 = 1 to z = 4. The left three panels are taken from 
IPorciani et al.l (120041 ) usin g the 2QZ d ata, whereas the right 
two panels are those from lShen et alJ (,2007) with the SDSS 
data. The marked redshift in each panel indicates the me- 
dian redshift of the measured quasars. For panels left-to- 
right, the absolute magnitude range defining each quasar 
sample is: the 2QZ survey Mb, £ [-25.32,-21.72], Mb, G 
[-25.97, -22.80], and Mb-, € [-26.44, -23.37], and the SDSS 
survejjfl Mi < —26.1 and Mi < —26.7 (all magnitudes have 
units of 5 log /i7o). 

The default model produces a very good fit to the ob- 
served quasar clustering at all redshifts considered. At 2 > 3 
there is a small under-prediction of the clustering amplitude, 
however the abundance of model quasars here is extremely 
low and the correlation function somewhat noisy. This over- 
all success gives us confidence that the use of the rriBH — cr 
relation and abundance mapping technique to construct our 
model is actually placing quasars in the correct halos, at 
least in a statistical sense. At 2 < 2 where the clustering is 



well measured, an upturn is seen at small scales 



IMpc). 



Hints of such excess pair s have been found i n the SDSS LRG- 
QSO cross c orrelation of Padmanabhan et al. (2008), as we ll 
as work bv iHennawi et all l|2006l ) and iMvers et al.1 (|2008l ). 
We leave a more detailed analysis of this result to future 
work. 



* Note that IShen et al. do not state the absolute magni- 

tude range that defines the quasars in their two high redshift bins. 
Hence, for simplicity we select model quasars brighter than the 
faint absolute magnitude corresponding to i = 20.2 (the SDSS 
quasar apparent magnitude limit) at the median redshift of each 
sample (M. Strauss, priv. comm.). 




Figure 9. Model predictions for the luminosity dependence of 
quasar bias, for quasar samples defined by a faint magnitude limit. 
Results at five redshifts are shown, from 2 = 0.5 to 2 = 4.0. 
Quasars at 2 ^ 1 show little or no luminosity dependent bias. At 
higher redshifts, a strong luminosity dependence in the bias is 
predicted for brighter quasars relative to faint. 



3.8 Luminosity dependent quasar bias at high 
redshift 

For most redshift bins our statistics are far superior to 
that which can be measured observationally. To test the 
model further we look for signatures of evolution in the 
clustering amplitude as a function of luminosity, captured 
here through the quasar bias. Such evolution, or lack 
there of, is a more detailed probe of how quasars oc- 
cupy halos, an d is currently the focus of much observa- 
tional scrutiny (iPorciani fc Norberg|200^ : lMvers et al ]|2007l : 



Ida Angela et al.ll2008h ." 



In figure P we show predict ions for the luminosity de- 
pendent bias fMo & White 1996) measured from our model, 
calculated using the Jenk ins ct al.. (,2001i ) mass function and 
a fit to the model quasar halo occupation fraction. The cal- 
culation is performed analytically to remove the noise of 
small number statistics from rare objects at the high mass 
end. Samples are defined with a faint limiting magnitude 
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(i.e. applying equation[7]to a halo mass cut), and we present 
results for five redshift bins ranging from z — 0.5 to z = 4. 

At z ^ 1 weak (or no) luminosity dependent bias is 
present in the model quasar population. At z > 1, how- 
ever, bright quasars show an increased bias with respect to 
faint quasars. This ranges from 6 ~ 3 to 6 ~ 7 between lumi- 
nosity extremes at 2; = 2, and fe ~ 5 to b ~ 15 at z > 2. Such 
increases are only marginally observable with current data, 
but will certainly be testable in future surveys as quasar 
numbers increase and the luminosity baseline widens. 

Both weak luminosity dependence at low redshift and 
significant luminosity dependence at high redshift constitute 
a firm prediction for the clustering of quasars in our model. 
These predictions arise as a consequence of the changing oc- 
cupation statistics of quasars in halos with time. At high 
redshift quasars are spread over a much wider range of halo 
masses relative to low redshift for a comparable (to L*) lu- 
minosity range (see figure I10|l . This produces the stronger 
clustering gradient seen in figure |9l Said another way, lumi- 
nosity dependent clustering evaporates at late times due to 
the narrowing of the range of halo masses that host quasars 
as the Universe ages. We will discuss this result further in 
Section mH 



4 DISCUSSION 

4.1 An (incomplete) overview of other popular 
quasar models 



Quasars have been modelled in a number o f 
in the pa st sev eral years JCiotti fc Ostrikei 



ways 



19971: 



Silk fc ReesI Il998l: iFabianI Il999l: lKa,uffmann fc Haehneltl 



2OOOI: iH achnch fc Kauffmann' I2O O0I: ICattaiied 



Ciotti_&p striker 
20031: iGranato et al 



2001 



Begelman fc Nath 



Pi Matteo et al.l 




Wvithe fc Loebl 
Kawata fc Gibsan[ 



Springel et al 



2001 



2002 



2005allbl: 



Hopkins et al.l 
Hopkins et al.l 



. .2005 
2005allb 



2005a 



2006a 



20061: iLidz et al.1 l2006l: iFontanot et al 



20061: IWvithe fc Padmanabhanll2006b: Malbon et all l2007l: 



Siiacki et al.ll2007l: iHopkins et al.ii2008a .bi: .Merloni fc Hein j 



20081 : iDi Matteo et all bOOSl ). The goal of such works has 
usually been to understand the cosmological evolution of 
quasars and their triggering mechanism and black hole gas 
accretion rates. This is typically achieved by matching the 
model output to observables like the quasar luminosity 
function, the evolving space density of bright quasars, and 
the TTiBH — o" relation. More recently, quasars (and more 
generally AGN) have been linked to the quenching of star 
formation in massive elliptical galaxies, and quasar models 
have adapted to refiect these new found appreciations. 

In fact, most modern AGN models are built upon 
the current belief that quasars ar e triggered by major 
merging events of ga s rich galaxies (|Di Matteo et al.ll2005l : 
iHopk ins et al]|2006al ). This is a reasonable assumption to 
make. Something significant must be happening to the gas 
in the galaxy to cause it to lose so much angular momen- 
tum, a necessary condition to drive gas into the central re- 
gions where the black hole resides. Without such angular 
momentum loss it is hard to imagine how the required near 
Eddington accretion rates can be achieved. 



From analytic arguments alone. [Silk fc Ree3(|l998l ) pos- 
tulate a critical black hole mass, determined by the sur- 
rounding halo properties, above which star formation has 
been suppressed due to an expanding quasar wind that 
sweeps the galaxy clean of its star forming gas. Black holes in 
this picture either form early in the collapsing proto-galaxy 
at above the critical mass, or form close the the critical 
mass and are maintained at this mass by the hierarchical 
growth of the system. This model provides an elegant ex- 
planation for many observed black hole - galaxy /halo corre- 
lations but does not predict other quasar properties such as 
their lumino s ities o r evolving space density. Regardless, the 
ISilk fc Reed (|l998t ) model has become somewhat of a seed 
from which a number of more detailed models have devel- 
oped. 

One popular e xtension of these ideas is disc ussed in a 
series of papers by IWvithe fc Loebl (IWvithe fc L ocb 200i, 
l2003f ). In their model, quasars and their subsequent rapid 
black hole growth are triggered from major mergers, as 
discussed above. Under the assumption that the local gas 
traps much of the quasar energy without radiating it away, 
and that the subsequent quasar luminosity is some fixed 
fraction of the Eddington luminosity, they derive a series 
of equations that allow them to predict quasar luminosity 
and black hole/host correlations at various redshifts. With a 
small number of free (but physically motivated) parameters 
their model is tuned to provide a good fit to the high red- 
shift quasar luminosity function, although it over-predicts 
the abundance o f bright low redshift quasars. The model of 
IWvithe fc Loebl is somewhat similar to ours but differs in 
one critical way. The abundance of quasars in their model 
is determined from halos who undergo rapid growth (i.e. 
major mergers). Our modelling makes no such assumption, 
but rather forces the correct quasar number by abundance 
matching to the actual quasar luminosity function. We dis- 
cuss the advantages of this below. 

To further explore this picture of black hole and galaxy 
growth a number of authors have turned to performing hy- 
drodynamic simulations of the complex merger and accre- 
tion processes themselves. It should be remembered, how- 
ever, that the physics of such processes cannot actually 
be resolved with current computing power, as the scales 
involved lie orders-of-magnitude below what is required. 
Never-the-less, a combination of the phenomenological and 
hydrodynamic methodologies allow for an increased level of 
detail and accuracy in the models that is unreachable by 
analytic methods alone. 

One exampl e of this is the work b y [Hopkins et al.l. In a 
series of papers (|Hopkins et al.ll2005al lbl, l2006allblld . l2007al lbl. 
;2008ai lbh. these authors explore the co-evolution of black 
holes and galaxies triggered by mergers using high resolu- 
tion hydrodynamic simulations which include many realis- 
tic physical processes. The most important aspect of their 
work is the inclusion of merger driven star formation and 
the redistribution of disk gas which drives both galactic 
bulge growth and, simultaneously, growth in the black hole 
through accretion. Once convolved with cosmological statis- 
ticfl e.g. the evolution of merger rates with time, their 
model produces results that match a large number of ob- 



iHopk ins et all simulate individual merger events, albeit a large 
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Figure 10. Halo virial mass predictions (solid lines) for quasars with measured luminosity and redshift, as defined by equation 1131 
Magnitudes are provided in both the 6j-band used by the 2QZ survey and i-band used by the SDSS (k-corrected to 2: = 2 - see 
equation [Toll. Dashed lines indicate the evolution of the host halo space density, as inferred from the observed quasar luminosity function. 
These contours show that the distribution of mass for quasar host halos narrows as the Universe ages. Horizontal dotted lines show 
the corresponding black hole mass, given by equation |6] The downturn of density contours relative to black hole mass at low redshifts 
indicate downsizing in the black hole population. Both upsizing and no-sizing are seen at higher redshifts, depending on the space density 
considered. 



servables: luminosity and mass functions, quasar lifetimes, 
Eddington ratios, and host galaxy properties. 

In a critical deviation fr om past models (including our 
work here), the simulations of lHopkins et al.l explicitly track 
the varying black hole accretion rate during the merger. This 
produces a range of predictions for how the quasar/AGN 
light curve should appear as a function of time, from first 
interaction to final merger remnant. Their simulations sug- 
gest a possible common origin for the many observed AGN 
types, a so c alled "unified" merge r driven model of AGN 
and galaxies (|Hopkins et al. I l2006al ). that depends primarily 
on the time during the merger that the system is observed. 
The beauty of their model is that it presents a well defined 
and testable picture of black hole and galaxy growth. The 
drawback to their work is that its detail and complexity 
sometimes cloud its interpretation (see Section [4.41 belowl . 



4.2 A tool for observers 

The quasar model presented in this paper provides a rela- 
tionship between quasar luminosity, redshift, dark matter 
halo mass, and black hole mass. The essence of this model 
is captured by equation [7] (Section I2.1.1|l . For convenience, 
we invert this equation to obtain 



log( 



Mvir 



1013 /i-i Mp 



(1.43 ±0.32) 



log (^j^ H{z] 



+ (0.72 ±0.11) log 



Lq/v 



1012 ft-^ Lq 



(13) 



where H{z) is defined by equation (2] Equation [13] can be 
used to take a quasar with observed bolometric luminosity 
Lq at redshift z and predict its halo virial mass, assuming 
values for parameters 7 (the relationship between virial and 
circular velocities) and rj (the Eddington luminosity fraction 
of the quasar). 

Equation [13] is plotted in figure 1101 From this figure 
alone the quasar host virial mass (either individual or aver- 
aged over a group) for a given quasar luminosity and red- 
shift may be simply read off (solid lines). We emphasise, fig- 
ure]!^ allows observers to determine statistically measured 
halo masses to complement their observations without the 
need for large quasar surveys or clustering measures. We 
discuss the uncertainty on such masses below. 

Also over-plotted in figure \TU\ are the corresponding 
mean space densities of quasar hosts (dashed lines), calcu- 
lated from the quasar luminosity functioijf] (Section l2.2p , and 
the black hole mass at fixed luminosity (horizontal dotted 
lines), as given by equation [5] 

The spacing of density contours relative to halo mass 
highlights a main result of this work, which is that the dis- 
tribution of the masses of dark matter halos hosting quasars 
narrows with decreasing redshift. This leads to behaviour 
such as luminosity dependent clustering at high redshift 



number of them, and not the evolution of structure in a cosmo- 
logical context. 



^ Note that the artificial bump at 2 = 3 is due to a change in 
the assumed fitting parameters of the quasar luminosity function 
(table [TJ and has no bearing on the results. 
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but not at low (see Section I3.8|l . and luminosity depen- 
dent quasar lifetimes at low redshift but not high (see Sec- 
tioning]). Also, as discussed Section [3. 5 1 the turnover of den- 
sity contours at z 2 relative to black hole mass demon- 
strates downsizing in the black hole population. At higher 
redshifts, rare massive black holes show no downsizing trend, 
whereas the more common L* quasars are predicted to be 
upsizing, especially at redshifts greater than z ~ 4. 

4.3 How well can halo mass be predicted? 

Before applying our model (particularly equation I13|l to an 
observation or set of observations it is prudent to understand 
the limits to which dark matter halo mass can be inferred 
given both the built in theoretical and observational uncer- 
tainties. 

Observational error is drawn from local measurements 
and has been propagated through each equation appropri- 
ately. Near the characteristic luminosity of the quasar pop- 
ulation, Lq/rj ~ IQ^^h^Q Lq, the mass error has magnitude 
0.32 dex in log units. One magnitude brighter or fainter than 
this increases the error to 0.34 dex, whereas two magnitudes 
translates to an error of 0.39 dex. While somewhat large at 
the extremes, this uncertainty is still sufficiently manageable 
that tight clustering constraints can be made and quasar 
host halo masses inferred from the model. 

Theoretically, our parameter 7 (describing the relation- 
ship between virial and circular velocities) can potentially 
take values of between 1.0 and 1.8, depending on halo mass 
and redshift. iPorciani et al.1 l|2004l ) argue that a value of 
unity is the most appropriate for high redshift quasars in 
the 2QZ survey, and hence is our choice. However, higher 
values may shift th e inferred ma ss down by up to 0.5 dex or 
more. According to lSeliakI (|2002l ). such a shift would primar- 
ily affect quasars hosted by lower mass halos than considered 
here, ~ IQ^^h'^MQ. 

Similarly, quasars are known to exhibit a range of 
Eddington values, somewhat in conflict with our single 77 
assumption (made to ke ep the model simple). However, 
iMcLure fc Dunlod (|2004l ) argue that r) — 0.5 is a reason- 
able mean value for the high redshift quasar population. If 
we instead assume an r] value of 1.0 we find a decrease in 
predicted halo mass by ~ 0.2 dex, while 77 = 0.1 results in 
an increase of predicted halo mass of ~ 0.5 dex. 

Across the entire quasar population we believe our pa- 
rameter choices are reasonable and justified (in the mean) by 
observation. However, individual quasars many not always 
be seen near their peak luminosity or have the virial-to- 
circular velocity ratios eissumed here. Under these circum- 
stances the mass inferred by our model will be incorrect. We 
emphasise that the quasar host halo masses predicted by our 
model are only accurate for the assumed 77 and 7 values, and 
within the measured observational error. Relative to these 
assumptions they should provide a valuable tool with which 
to probe the quasar population across a wide luminosity and 
redshift range. 

4.4 Benefits to keeping it simple 

So what advantages does our model provide over past works? 
First, and despite much circumstantial evidence (a lot of 



which is very convincing) , astronomers still do not know the 
actual conditions and caveats under which quasar triggering 
occurs. We may be wrong about the merger hypothesis, or 
the circumstances under which the triggering is otherwise 
(in)effective, or there may exist more than one mechanism to 
trigger quasars. By using the data to determine which halos 
host quasars we free ourselves from pre-deciding the quasar 
trigger and the (perhaps unappreciated) consequences this 
may bring. 

Second, it is useful to build a statistically accurate rep- 
resentation of the quasar population (or at least one repre- 
sentation), from which we can 'work backwards', so to speak. 
Once our model is constructed we can examine it with confi- 
dence knowing that it has been constrained to be correct. A 
statistically correct model of the quasar population also al- 
lows observers to explore non-physics related issues in their 
data and survey design, such as cosmic variance and system- 
atics. 

Third, our model is largely transparent in its cause and 
effect, which makes it easy for both theorists and observers 
to understand and apply. This is rather a critical point. One 
of the primary applications of any theoretical model is as 
tool to interpret the data in a physically meaningful way. 
Unless it is clear why a model behaves the way it does there 
is little insight to be gained from matching the data alone. Of 
course, it is through the marriage of techniques that are both 
simple (to build intuition and set direction) and detailed (to 
understand the actual physics) that progress is made. 



5 SUMMARY 

Under minimal assumptions we have demonstrated a (statis- 
tically) successful phenomenological method to occupy dark 
matter halos with quasars in a way consistent with many 
key observations out to z ~ 5. We summarise the primary 
results and specific predictions of the model: 

• We provide simple equations to predict host dark mat- 
ter halo mass (eauation ll3p and black hole mass (equation[8} 
for quasars of given luminosity and redshift. These equations 
are applicable to both single and group quasar observations. 

• We pro vide a new j o int fi t to the qu asar luminosity 
functions of ICroom et al.l l|2004l) (2QZ) and [Richards et al.1 
(2006) (SDSS) which is accurate across QA < z < 4.8 (Sec- 
tion [33] and figures OH]) . 

• High redshift quasars of a given luminosity live in less 
massive dark matter halos than quasars of the same lumi- 
nosity at low redshift. Another way to say this is that, at 
fixed halo mass, high redshift quasars are brighter than their 
low redshift cousins (Section [STT] and figure [l]). 

• Our model predicts luminosity dependent quasar life- 
times of tq ~ 10^~*^7Q^yr, but which may be as short 
as 10^~^ yr for quasars brighter than L* and z 2. At 
z ^2 this bright trend is less pronounced and even reverses 
(Section 13.41 and figure [S]). This occurs because low red- 
shift bright quasars are rare amongst the fairly abundant 

mass halos that host them, whereas high redshift 
bright quasars are frequent among the rare 10^"^ mass 
halos they occupy (figure [TOj. 

• Our active black hole mass function is consistent with 
that observed but predicts a significant excess (factor of 5) 
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of fainter quasars (logmBH ~ 8.5) at z ~ 2 (Section 13.51 and 
figure |6]) . 

• "Downsizing" naturally arises in our model; at fixed 
space density black hole mass decreases significantly at red- 
shifts less than 2, while at higher redshifts this downsizing 
trend disappears for low space density objects and even re- 
verses (i.e. "upsizing") for higher space density objects (fig- 
ures |6] and [lOl) . 

• We predict evolution in the amplitude of the ttibh — 
A^haio relation with time, with black holes of a given mass 
increasingly hosted by less massive halos at earlier times. 
The amplitude changes by about a factor of 20 between z — 
5 and z = (Section 13.61 and figure [T]). 

• Our model quasars have the correct clustering proper- 
ties when compared to observations out to z ~ 4 (clustering 
was not a constraint on the model). However, our clustering 
amplitude may be slightly low for quasars at 2 ^ 3 (Sec- 
tion [37] and figure [8]). 

• Our model places quasars in halos in such a way that 
very little (or no) luminosity dependent clustering exists at 
z ^ 1 for all magnitudes. Ifowever, we predict strong lumi- 
nosity dependent clustering at higher redshifts for luminous 
quasars when compared with the L* population (Section l3.8l 
and figure[9]). This behaviour results from the narrowing dis- 
tribution of halo masses that quasars occupy as the Universe 
ages (figure [TUl) . 

Quasars appear set to remain a valuable probe of galaxy 
formation out to high redshift due to the distances they can 
cleanly be measured. They are furthermore a incredibly in- 
teresting population of objects in their own right, encapsu- 
lating significant amounts of fundamental physics and broad 
phenomenology still yet to be understood. The ability to 
produce statistically accurate models of this unique popula- 
tion will be essential to interpreting their future observation. 
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